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Abstract. We study the effects of diffusing opinions on the Deffuant et al. model for continuous opinion 
dynamics. Individuals are given the opportunity to change their opinion, with a given probability, to a 
^ randomly selected opinion inside an interval centered around the present opinion. We show that diffusion 

induces an order-disorder transition. In the disordered state the opinion distribution tends to be uniform, 



while for the ordered state a set of well defined opinion clusters are formed, although with some opinion 



spread inside them. If the diffusion jumps are not large, clusters coalesce, so that weak diffusion favors 
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opinion consensus. A master equation for the process described above is presented. We find that the master 
equation and the Monte-Carlo simulations do not always agree due to finite-size induced fluctuations. Using 
a linear stability analysis we can derive approximate conditions for the transition between opinion clusters 
and the disordered state. The linear stability analysis is compared with Monte Carlo simulations. Novel 



^1- 

interesting phenomena are analyzed. 

O 

> 

PACS. 8 9.65.-S Social and economic systems., 05. 40.-a Fluctuation phenomena, random processes, noise, 
3 and Brownian motion. 

1 Introduction ity of individuals adopt a similar opinion (consensus) or 

a number of opinion groups (clusters) arise. In many sit- 
in a community, opinions evolve due to affinities and con- uations, the dynamics of this complex collective behavior 
diets between mutually interacting individuals. These in- goes beyond specific individual attributes and seems to 
tcractions lead to collective states, where either a major- be well characterized by quantities like statistical distri- 
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butions and averages. This explains why, in the last years, model in which we allow individuals to change their opin- 

the understanding of opinion formation using tools and ion, with a given probability, to a randomly chosen value 

techniques borrowed from nonlinear and statistical physics in the whole opinion space. One possible interpretation of 

has become a topic of interest for physicists. In particular, this modification (which, from the technical point of view, 

several models have been developed to reproduce the ba- can be understood as some sort of noise acting upon the 

sic elements that drive the processes of opinion evolution dynamics) is that individuals keep at all times a basal 

PP. These models can be classified in two broad groups: on opinion to which they return from time to time, no mat- 

the one hand, discrete, Ising-type, models where opinions ter what the environment tells them to do, modelling in 

can only adopt a finite set of integer values 2,3; on the a crude way the intrinsic free-will of human decisions. We 

other, continuous opinion models where opinions can vary have shown [13] that this modification of the dynamical 

continuously in a finite interval 4,5,6,7,8,9,10 . rules can induce new and interesting phenomena, such as 



noise-induced bistability, and that the model displays in- 
teresting and non-trivial finite-size effects that need to be 
considered in some detail. 



In the context of continuous opinion dynamics a model 
introduced by Deffuant and collaborators has received much 

attention recently 6J. In this model individuals meet in 

, . . j 1 1 , In this paper, we continue our analysis of the uncer- 

random pairwise encounters and then converge to a com- r r 1 J 

i .j- ,i . , . . . • , •, tainty factors present in human decisions by considering a 

mon opinion only it their respective opinions prior to the jo 

,-rr 1,1 • i • i r modification of the model in which the opinion of an indi- 

encountcr diner less than some given amount, m a kind ot 1 

, , , ni i h i i n .m »n , vidual can iump to a new value chosen randomly inside an 

bounded confidence mechanism [111112] . After some tran- j r j 

, ,. in . , , , n , . , . , .,, interval of finite width centered around the current opin- 

sient evolution, this leads to hnal states m which either 

r . , n , ,. ,., . n ion. In other words, in addition to the bounded confidence 

lull consensus is reached or the population splits m a h- 

, r i i n , n*j**j i mechanism, the opinions execute a random-walk motion 

mte number of clusters such that all individuals m one 

, , tt i_ v ii „ or diffusion in opinion space. When the typical size of the 

cluster share the same opinion. However, we believe that r 

, , . . ,, r , , , , r . iumps is sufficiently large we expect this model to behave 

such dynamics misses the fact that the behavior of an m- J r jo r 

n , , . a as the previous free-will model discussed in 1 13 , but when 

dividual m a society does not depend only on the influence 1 — ' 

r i .,i ■ -1 mi jj-i- i it is small a distinct behavior will be obtained, 

of people with similar opinions. Ihere are many additional 

factors, both personal and external, which also cause opin- We analyze which aspects of the Deffuant et ad. model 
ion changes. These factors may introduce some degree of are robust against the introduction of this diffusion pro- 
discrepancy within otherwise well defined opinion clusters cess. The analysis, based upon Monte Carlo simulations 
and should be modeled somehow. In a recent paper [13] as well as on numerical integrations of the corresponding 
we analyzed the effects induced by a modification of the master equation, reveals indeed a dependence on the size 
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of the random jumps. For small jumps, we have found formation of opinion clusters; summary and conclusions 
that the center of mass of each opinion cluster executes a are presented in section [4j 
random walk with an effective diffusion coefficient equal 
to the one corresponding to a single individual divided 
by the number of individuals in the cluster. The interplay 
between interactions and random walk induces a coarsen- 
ing process that leads to the formation of a single large 
cluster at very long times. When large jumps are allowed, 
the clusters tend to stay in relatively well defined posi- 
tions. We have observed, in a way similar to the model 
in [15] . regions of bistability in which fluctuations can 
take the system from one state to another and back. We 
will show that Monte Carlo simulations do not necessarily 
agree with the results of the master equation because of 
the inherent finite-size-induced fluctuations. Finally, our 
studies reveal that there is a region in parameter space in 
which the system becomes disorganized and cluster forma- 

tion does not occur. We provide a linear stability analysis x n+i ~ x n+i = n > (v 



2 A continuous-opinion dynamics model with 
bounded opinion jumps 

The original version of the Deffuant ct ad. model [6] con- 
siders a population with N individuals. The opinion x l n 
on a given topic that individual i has at time-step n is a 
real variable in the interval [0, 1]. One also assumes that 
the initial values x l for i — 1, . . . , N are randomly dis- 
tributed in this interval. A bounded confidence mechanism 
is introduced to reflect that individuals interact, discuss, 
and modify their opinions: at time-step n two individuals, 
say i and j, are randomly chosen; if their opinions satisfy 
\ x n ~ x n\ < e > they converge to the common value: 



that qualitatively reproduces this order-disorder transi- 
tion. Analytical results are compared with Monte Carlo 
simulations using the so-called cluster coefficient which 
aims to characterize the existence of clusters. 



otherwise they remain unchanged. Whether the opinions 
have been modified or not, time increases n — > n+ 1. It is 
customary to introduce the time variable t = nAt, where 
At = 1/N, measuring the number of opinion updates per 
individual, or number of Monte-Carlo steps (MCS). As a 
This paper is organized as follows: the model and main consequence of the iteration of this dynamical rule, the 
results coming from Monte Carlo simulations and numer- system reaches a static final configuration, which depends 
ical integration of the master equation are presented in on the confidence parameter e taking values between 
section[2j in section[3]we use the so-called group coefficient and 1. Starting from uniformly distributed random values 
to characterize the order-disorder transition that appears for the initial opinions, the typical realization is that for 
in the model and we apply a linear stability analysis to e > 0.5 the system evolves to a state of consensus where 
derive approximately the critical parameter values for the all individuals share the same opinion and that, decreasing 



4 M. Pineda, R. Toral, E. Hernandez-Garcfa: Diffusing opinions in bounded confidence processes 

e, the population splits into opinion clusters separated by more realistic, adsorbing boundary conditions, in which 

distances larger than e [STTT4"]. opinions that try to go away towards the left or towards 

the right of the interval [0, 1] are set to or 1, respectively. 

The new ingredient we add to the dynamics is that 

Nevertheless, for mathematical simplicity periodic bound- 
individuals can perform random jumps in their opinions. 

ary conditions will be considered in some particular cases 

More specifically the dynamical rules are modified as tol- 
as properly mentioned. 

lows: at time step n the dynamical rule Eq. (JlJ applies 
with probability 1 — to; otherwise a randomly chosen in- 

2.1 Monte Carlo simulations 

dividual i changes the opinion to a new value x l n+l ran- 
domly chosen from the interval (x l n — "/,x l n + 7). The pa- We present in this subsection the main results obtained 
rameters 7 € [0, 1] and to € [0, 1] determine, respectively, from Monte Carlo simulations for a system with a finite 
the width and frequency of random jumps. When 7 is population of N individuals and adsorbing boundary con- 
large enough, we expect that the behavior of this model ditions. As in previous studies, we assume that the ini- 
will approach the one of our previous free-will model |13j , tial condition represents a uniform distribution in opinion 
in which noisy jumps occurred to random locations in space interval [0,1]. In the original Deffuant et al. model 
opinion space. Since the variance of the opinion jumps (to = 7 = 0), Monte Carlo simulations show that for e > 
is a 2 — 7 2 /3, each individual would behave as a random and starting with homogeneous initial conditions the sys- 
walker with diffusion coefficient D = to7 2 /3 if interac- tern either reaches a final state of perfect consensus or 
tions and boundaries were absent. Note, however, that it splits into a finite number of clusters such that all indi- 
is possible from the model rules that opinions leave the viduals in one cluster have exactly the same opinion [6]. 
bounded opinion space [0, 1]. In order to assure that opin- This picture changes when diffusion is introduced in the 
ions stay inside this interval, we need to implement proper way described above. First, for to > and 7 > 0, clusters 
boundary conditions to the dynamics. A simple and math- still are formed but individuals in a single cluster do not 
ematically convenient choice is to use periodic boundary have exactly the same opinion but there is some disper- 
conditions, so that the opinion interval is considered to be sion among them. Second, when 7 is small the center of 
wrapped on a circle. In this case, opinions that go away a mass of each cluster performs a random walk through the 
certain distance to the left of the extreme opinion or to opinion space. As a consequence, opinion clusters start to 
the right of the extreme opinion 1 are injected throughout collide and merge in a coarsening process that leads li- 
the opposite extreme by a similar amount. In real situa- nally to a single surviving large group at very long times, 
tions it is unlikely for individuals with extreme opinions To visualize these interesting behaviors, Fig.[T]shows time 
to change their opinions so drastically. Thus, we will use series of the opinion space for a value of e such that only 
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one big cluster is formed. Dispersion around the center of 
mass remains approximately independent of N, but it is 
seen that the center of mass X cm of the cluster behaves as 
a random walker in opinion space, with smaller displace- 
ment for larger N. We can understand this motion by 
noting that the binary opinion interactions do not change 
the position of the center of mass, which then moves only 
because of the diffusive motions of the opinions. We can 
thus write for the center of mass X cm at time t + 1 (i.e., 
after the number of steps has been increased from n to 
n + N) as its position at time t plus the motion induced 
by the individual jumps occurring in that time interval 
(let us say that there have been N\ of them): 

N Nt 
X cm (t+1) EE jj^X^t+l) = ^cmW + ^X)^*)- (2) 

i=l i=l 

£*(t) are independent random jumps of zero mean and 
variance a 2 = 7 2 /3. Reordering, and taking squares and 
mean value one obtains 

(ix cm (t+i)-x cm (t)} 2 ) = /f^E^wl ) 

mN 2 mj 2 

= Jp a = iVT' (3) 
where we have used that the expected value of Ni is mN. 
This implies that the cluster's center of mass experiences 
a random walk with an effective diffusion coefficient D cm 
equal to the single-individual's opinion one divided by the 
number of individuals in the cluster D cm = D/N . Fig. [T] 
illustrates cluster random walks for three different popula- 
tion sizes. For isolated clusters (i.e. separated by distances 
larger than e) the arguments above are only valid when 
boundary effects are unimportant, which means that clus- 
ters remain sufficiently far from the interval extremes (at 



a) 




0123-1 1234 1234 
t (10 4 MCS) t (10 4 MCS) t (10 4 MCS) 



Fig. 1. Time series of the opinion distributions at three values 
of TV. m = 0.1, 7 = 0.1, and e = 0.28. (a) At N = 100 the 
cluster moves around the whole opinion space, (b) The cluster 
at N = 1000 moves less, (c) At N = 10000 the cluster remains 
close to 0.5. Simulations are performed with homogeneous ini- 
tial conditions and adsorbing boundary conditions. Although 
the simulations were done with 100, 1000 and 10000 agents, 
here we plot in all panels only 100 of them, and at intervals of 
100 MCS, to avoid saturation of the plots. 

distances larger than 7) so that jumping particles can not 
reach the boundaries. For larger 7 or e particles will feel 
the boundaries more easily and we expect cluster mobility 
to be much reduced by boundary effects. 

From Monte Carlo simulations, it is also seen that clus- 
ter formation always occurs for small values of e, 7, and m. 
But, in contrast with the diffusionless model, coarsening is 
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t (10 5 MCS) 



Fig. 2. Time series of the opinions with e = 0.05, m = 0.01, 
7 = 0.04, and N = 1000. To avoid saturation of the plot, only 
100 agents are shown, and at intervals of 450 MCS. The clusters 
perform random walks and successive merging of clusters oc- 
curs after collisions. Finally, at very long times, a state contain- 
ing a single cluster will dominate. Simulations are performed 
with homogeneous initial conditions and adsorbing boundary 
conditions. 

observed for small 7. Figure [2] shows successive merging of 
clusters, occurring after the collisions which arise because 
of the diffusive wandering of clusters (by collision we mean 
that clusters become closer than e, so that they interact). 
We expect that in this regime, at very long times, a state 
containing a single cluster would be the final regime. 

The behavior is different at large values of 7: in ad- 



which tend to stay at relatively well defined positions, we 
see that there is no tendency to reducing the number of 
clusters (See Figs. [3^, and c), so that a pattern of opinions 
is established. The pattern of clusters is approximately pe- 
riodic when £ < 1. This is similar to the behavior of the 
model in [T3]. In fact, for 7 of the order of system size the 
jumps are of global size so that there should be no differ- 
ence between both models. As in [13], a remarkable fact 
occurs in which one can find regions of bistability: inside 
these regions, the inherent fluctuations arising from the fi- 
nite number of particles take the system from one state to 
another and back. The sort of transitions between steady 
states is observed for instance, in Monte Carlo simulations 
at e = 0.316. Figure [3]d shows multiple jumps between a 
state of two big opinion cluster and another state of a big 
cluster with two smaller ones near the edges of the opinion 
interval. 

2.2 Master equation approach 

In this subsection we analyze the master equation descrip- 
tion of the process introduced above. Following standard 
arguments (see for example [13) ) one finds the master 
equation for the probability density function P{x, t) of an 
individual opinion x at time t 
dP(x,t) 



Of 



\x~x 2 \<t/2 



(1 - m)[4 



-2P(x,t) / dx 2 P{x 2 ,t)] 

J \x— x 2 \<e 

+m[G(x,t) - P(x,t)] . 



dx 2 P{2x - x 2 ,t)P(x 2 ,t) 



(4) 



The term proportional to (1 — m) is the one coming 



dition to the expected reduced wandering of the clusters, from the original rules of the Deffuant et al. model 0Q3 
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boundary conditions: 
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Fig. 3. Time series of the opinion distributions at three values 
of e and m = 0.01, 7 = 0.4, and N = 1000. To avoid saturation 
of the plot, only 100 agents are shown, and at intervals of 1000 
MCS. (a) At e — 0.28 the system appears polarized in two 
opinion clusters, (c) A single major cluster with two lateral 
minor clusters is observed at e = 0.35. (b) At e = 0.316 the 
systems fluctuates between these two states. Simulations are 
performed with homogeneous initial conditions and adsorbing 
boundary conditions. 



whereas the one proportional to m describes the random 
jumps. The function G(x, t) is, for the case of adsorbing 



if x < 7, 



if 7 < x < 1 — 7, 



5(x-l)fl 7 dx'=±%±^P(x',t) 
+ SLrW P ( x '> t '>> 



if x > 1 — 7. 

(5) 



For small values of 7 the boundary effects become less 
important, and the second case in |5| applies in the ma- 
jority of cases. In addition, for 7 small enough, this term 
can be approximated as 



x+1 dx' * , , . 7 2 d 2 P(x,t) 



(6) 



so that the diffusion term, i.e. the term proportional to 
to in the right-hand side of Eq. Q becomes of the form 
Y 9 g^'^ , where one finds again the diffusion coefficient 

D = 7777 2 /3. 

We have solved numerically the master equation Q for 
the distribution P(x, t) starting from an initial condition 
representing an uniform distribution in opinion space, i.e. 
P(x,t = 0) = 1 for x € [0, 1] and P(x,t = 0) = oth- 
erwise. For m = it is well known that the distribution 
Poo{x) = lim t _>. 00 P(x, t) is a sum of delta- functions lo- 
cated at particular points [5,14 . However, this is not the 
case for the full Eq. Q with m > 0: the final distribu- 
tions are no longer made up of delta-functions and their 
final shape strongly depends on the particular values of 7, 
e, and m. We have found that for small values of 7 only 
one opinion cluster, centered around x = 0.5, remains for 
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almost any value of the parameter e (see Fig.[4p,). We in- 
terpret this as a consequence of the dynamics discussed in 
section |2.1[ where we showed that although several clus- 
ters were initially formed, collisions reduced their number 
and a single one was expected to survive at long times. 
We stress that the single cluster in Fig. [4^, is only ob- 
tained at very long times. At shorter times several clusters 
are present in the solution of the master equation start- 
ing from a flat initial condition (see for example Fig. [5]). 
The width of the cluster in this master equation descrip- 
tion should be related to the dispersion seen in the Monte 
Carlo simulations (see Fig. [I]) at a given time, and not 
to the amplitude of the diffusive cluster wandering, since 
the master equation description is expected to be accu- 
rate as N — > oo, a limit in which cluster Brownian motion 



becomes frozen [see (Eq.pl) and section 2.3 



For large values of 7 the situation is rather different: as 
in the model with free-will |13] , for small e several opinion 
clusters form, which do not coalesce. Asymptotic states as 
a function of e are displayed in Fig. |3]d. Different cluster 
bifurcations, as in the original Deffuant model [5] and in 
the case of free-will [T2] are seen. 

An important feature is that, for all values of 7 con- 
sidered, the clusters become less defined below a critical 
value of e (for large 7 one can always observe however the 
two large clusters at x = 0, 1 arising from the adsorbing 
boundary conditions) . This point will be further addressed 
in section [3l 




7=0.1 



7=0.4 



Fig. 4. Plot of the asymptotic probability density P(x, t — > 
00), coded in logarithmic grey levels, as a function of e, ob- 
tained after a numerical integration of Eq. (|4| starting with a 
flat distribution. The top panel is for 7 = 0.1, and the bot- 
tom panel for 7 = 0.4. In both cases we used m = 0.1 and an 
integration time t = 5 x 10 4 . 

2.3 Master equation description versus Monte Carlo 
realizations 

Monte Carlo simulations with a finite number N of in- 
dividuals are not always faithfully represented by the so- 
lutions of the master equation, which implicitly assumes 
that N — > 00, in addition to considering correlations be- 
tween individuals only in an approximate way. In Fig. [5] 
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Fig. 5. Probability distribution function P(x,t), for interme- 
diate time steps, from Monte-Carlo (MC) simulations (his- 
tograms binned with bin size Ax — 5 x 10 -4 ) for two different 
system sizes N = 10 3 and N = 10 5 , and the master equation 
(ME) integrations of Eq. Q starting with a flat distribution 
with e = 0.28 and m = 0.1. The Monte Carlo distributions are 
an average over 10 5 realizations for N = 10 3 and 10 4 realiza- 
tions for N — 10 . The top panels present results with 7 = 0.1 
and the bottom panels with 7 = 0.4. 



we plot the time evolution of the probability coming from 
Monte Carlo simulations (an average of the binned opin- 
ion distribution over a large number of realizations is per- 
formed, see caption) and the results from the master equa- 
tion in the case m — 0.1 and e = 0.28, which is close to a 
bifurcation from one to two big clusters in the diffusionless 
model (see [5] for details). The top panels correspond to 
7 = 0.1 and the bottom panels to 7 = 0.4. It can be seen 
that, although the Monte Carlo simulation and the mas- 
ter equation agree initially very well, they start to deviate 
after a time that depends on the number of individuals 
N: the larger N, the longer the time for which the Monte 
Carlo simulations are faithfully described by the master 
equation. 

In the case 7 = 0.1, while the numerical solution of the 
master equation and Monte Carlo simulations both reach 
at long times a distribution with a single large maximum 
(large cluster) at x = 0.5, in the Monte Carlo simulations 
the width of this steady distribution is larger for small 
particle number N . One can understand this by noticing 
that the Monte Carlo result is in fact the average over a 
large number of realizations, and each one of them consists 
on a single cluster whose location fluctuates widely (see 
Fig. [I]) . The fluctuations in the location of the center of 
mass of the cluster are reduced for increasing N (see sec- 
tion 



2.1 ), so that only the natural width of the cluster will 



show up for increasing JV, the range in which the mas- 
ter equation description is expected to be accurate. For 
7 = 0.4, we recover the same type of behavior already re- 
ported in our previous work concerning unbounded jumps 
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in opinion space [13]. One can see in the bottom panel of 3 Order-disorder transitions 

figure [5] that while the numerical solution of the master _ 

As Fig. Ml shows, for e smaller than a critical value (which 

equation tends at long times to a unimodal steady-state 

depends on m and 7) the probability distribution becomes 

distribution with a single large peak (surrounded by two 

blurred such that the maxima of the distributions are not 

small ones), the Monte-Carlo simulations end up at long 

evident, implying the inhibition of cluster formation and 

times in a bimodal one with two peaks. At intermediate 

the establishment of a more homogeneous opinion dis- 

times both types of solution agree, the larger the value 

tribution (although there is always some inhomogeneity 

of N the longer the time of agreement. A similar discrep- 

close to the boundaries, specially for large 7). A similar 

ancy between the Monte-Carlo and the master equation 

effect can be observed with Monte-Carlo simulations un- 

results also appeared in our previous study concerning un- 
der adsorbing boundary conditions and can be described 

bounded iumps[T3"]: starting from a uniform initial condi- 

in terms of an order-disorder transition: order identified 

tion, a unimodal steady state is reached by the master 

with the state with well defined opinion clusters and dis- 

equation but a bimodal distribution is the one reached in- 
order identified with the state without clusters. 

stead in Monte-Carlo simulations. The point is that both 

To identify in a more quantitative way this order-disorder 

distributions (the unimodal and the bimodal) are station- 

transition, we use the so-called cluster coefficient Gm [131 

ary solutions of the master equation, but the unimodal 

[To] . Its definition starts by first dividing the opinion space 

solution is metastable: a perturbation would take the sys- 

[0, 1] in M equal boxes and counting the number of indi- 

tem out of this solution towards the bimodal distribution. 

viduals li which, at time step n, have their opinion in the 

The perturbation needs to break the left-right symmetry 

box [{i—l)/M,i/M]. The value of M must not be so large 

(or x — > 1 — x) of the problem which is present for a 

that particles are artificially considered to be part of a 

uniform initial condition. In the case of the Monte-Carlo 

single cluster, nor so small that statistical errors are large 

simulations the perturbation is induced by the unavoid- 



able finite-size fluctuations, thus explaining the discrepan- ^ M ; j 



within one box. We choose M — 100. One next defines an 
entropy Sm = — X)f=i n m N> an< ^ ^ ne cms ter coefficient 



cies. In the case of the numerical integration of the master __ 

G M = M- 1 (e s ™), (7) 



equation, this kind of perturbation appears if one is not _ , ■ 1 / s 



careful enough and introduces, for example, round-off mi- 

where the over-bar denotes a temporal average in steady 

merical errors that do not respect the above-mentioned 

conditions and (•) indicates an average over different real- 
symmetry. 

izations of the dynamics. Note that 1/M < Gm < 1. Large 
values, Gm ~ 1, indicate that the opinions are evenly dis- 
tributed along the full opinion space (a situation identi- 
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fied with disorder), while small values of Gm indicate that 
opinions peak around a finite set of major opinion clusters 
(a situation identified with order). 





Fig. 6. Opinion cluster coefficient Gm versus m for 7 = 0.1 (a) 
and 7 = 0.4 (b) obtained from Monte Carlo simulations with 
N = 10 4 and adsorbing boundary conditions. The values of 
the confidence parameter are e = 0.05, 0.2, 0.4, and 0.6, from 
left to right (dots; the solid line is a guide to the eye). We will 
define the transition from order to disorder as the location m c 
of the absolute maximum value of Gm- 

In Fig. [6|a) we plot Gm as a function of m for 7 = 0.1 
and different values of e as obtained in the Monte Carlo 
simulations using adsorbing boundary conditions. For this 



value of 7, Gm reaches an absolute maximum value close 
to 1 (corresponding to a completely unstructured state) 
and then decreases monotonically. The adsorption by the 
extreme opinion values and 1 prevents the formation of a 
fully homogeneous state (which will give Gm = 1) as two 
opinion clusters are formed at the extrema of the opinion 
space, whereas the rest of the opinion space remains more 
or less homogeneously populated. Fig. [6^6) shows a similar 
behavior in the case 7 = 0.4, but in this case the number 
of individuals whose opinion is adsorbed by the extremes 
is larger. Therefore, Gm saturates farther away from its 
maximum possible value 1. We will define the transition 
from order to disorder as the location m c of the absolute 
maximum value of Gm- In the next subsection we will 
explain this transition via a simple linear stability analysis 
that turns out to be very accurate, in particular, for small 
values of 7. 

3.1 A linear stability analysis 

Although the transition to cluster formation is a nonlinear 
process, one can still derive approximate analytical condi- 
tions for the existence of cluster formation as a function 
of the control parameters by performing a linear stabil- 
ity analysis of the unstructured solution of Eq. Q. This 
is greatly simplified if one neglects the influence of the 
boundaries and assumes that there are periodic boundary 
conditions at the ends of the interval [0, 1]. This would be 
a reasonable approximation for describing the distribution 
far from the boundaries. We expect this approximation to 
be valid for not too large 7 or e since, as seen for exam- 



Y=0.1 
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pie in Fig. |4j there is not much structure near the edges 
of opinion space in those cases. Under periodic boundary 
conditions, the homogeneous configuration Ph(x) — 1 is 
the unstructured steady solution of the master equation 
which results from Eqs. (4|5) without the contribution of 
the boundary terms: 



7=0.4 




dP(x,t) 
dt 



= (1 - m)[4 / dx 2 P(2x - x 2 , t)P(x 2 ,t) 

-2P(x,t) f 

J \x — x 2 \ <e 

(8) 



in 



1 


2 3 4 5 




qt 


:(b) 









0.2 0.4 0.6 0.8 

e 




' dx 2 P(x 2 ,t)} 

J \x— x 2 \ <e 

X+7 ^P^t)-P(x,t) 

x 7 J Fig. 7. (a) Growth rate, equation (|9j), of a perturbation to the 

To analyze the stability of the homogeneous solution homogen eous state as a function of qe for e = 0.23, 7 = 0.1, and 

m = 0.0, 0.4, and 0.8, from top to bottom, (b) Phase diagram 



P(x) = 1 we write P(x,t) = 1 + A q e lqx+Xqt , where q is 



the wave number of the perturbation, X q its growth rate on the plane (e, m) as predicted by the linear stability anal- 
and A q the amplitude. Introducing this ansatz in Eq. (|8| ysis with 7 = 0.1 (solid line) is compared with Monte Carlo 



we find the dispersion relation 

4sin((je/2) sin^e) 



Xq = 4e(l - m) 



+m 



sin(<n) 
97 



1 



(9) 



We note that, for small 7 the second term, propor- 
tional to m, becomes —m{ r )q) 2 /Q as corresponding to the 
expected diffusive behavior. The top panels of Fig. [7] show 
this dispersion relation as a function of qe for several val- 



simulations using adsorbing boundary conditions and N = 10 4 
(open dots). Disordered states are above the lines and ordered 
(clustered) ones below, (c) Growth rate as a function of qe for 
e = 0.23, 7 = 0.4, and m = 0.0, 0.1, 0.15 and 0.2, from top 
to bottom, (d) Phase diagram on the plane (e, ?n) as predicted 
by the linear stability analysis with 7 = 0.4 (solid line) is com- 
pared with Monte Carlo simulations using adsorbing bound- 
ary conditions and N = 10 4 (open dots) . Disordered states are 
above the lines and ordered (clustered) ones below. 



ues of m with 7 = 0.1 and 0.4. If there exists a wavenum- 
ber q for which X q > 0, then the unstructured uniform ical expressions can be obtained expanding X q in powers 



state is unstable and cluster formation will be possible, of q 
It is not possible to find closed expressions for the wave- 
length q max giving the maximum growth rate Xq max , or 
the critical values m c and e c defining the regions where 
the homogeneous state is stable for fixed 7, but all of this 



X q = 



mrf (1 - m) 3 
6 3 6 



1 

120 



777,7 



7(1 - m)e 5 



q 4 + 0(q 6 ). (10) 



In the case in which the q 4 term remains negative 
can be readily obtained numerically. Approximate analyt- (which occurs if the 7777 term remains smaller than the 
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one containing e), the change of sign of the q 2 term iden- dieted by the linear stability analysis. We see in the figure 

,.,2 ll/3 



tifies e r = 



as the value of e below which an that there is a very good agreement between theory and 



2(l-m) 

unstructured configuration remains stable. Alternatively, simulations for 7 = 0.1, and a worse correspondence, al- 

for fixed e we find m c = 2e ^+j 2 > ^ ne critical value above though still qualitatively correct, for 7 = 0.4. Thus we can 

which clusters will disappear. Within this approximation say that the arrest of cluster formation arises because the 

the expression for the fastest growing mode near the order- random jumps stabilize the homogeneous opinion distri- 

disorder transition is: bution. 



(e-e c ) . (11) 



4 Summary and conclusions 

We stress that all these expressions following Eq. ( 10 1 

are only valid for the case in which the appearance of We have studied the continuous opinion model by Def- 

positivc values of X q when varying a parameter occurs fuant efc &1. when one adds diffusion to the dynamics. More 

first at values of q close to zero, corresponding to a long- precisely, we have modified the evolution rules by includ- 

wavelength instability. As seen in Fig. [7] this is not always ing the probability that individuals change their opinion 

the case for large values of 7, for which the transitions to another value randomly chosen inside an interval cen- 

have to be obtained numerically. In fact, we note that tered around the current opinion. Our aim is to include the 

the change in behavior observed in the simulations from effects of the unavoidable elements of randomness always 

a tendency to coarsening to the establishment of a robust present in human decisions. When the typical size of the 

periodic pattern, occurring for increasing 7, seems to be random jumps is large, this model approaches a previous 

correlated with the change in the character of the linear minimalistic model introduced by us [13]. 
instability, from long-wavelength to finite wavelength, this The final collective states depend on the system size N, 

last one being the one occurring in the free-will model in the confidence parameter e as well as on the typical size 

|13j . A more detailed discussion of the relationship be- of the random jumps 7 and their probability m. While 

tween the linear dynamics and the nonlinear long-time our numerical results consider the more natural adsorbing 

states is however beyond the scope of this paper. Compar- boundary conditions in which opinions beyond the limits 

ison of the predicted order-disorder transition with Monte of the allowed interval are set to the extreme values or 

Carlo simulations is performed in the bottom panels of 1, some analytical calculations are carried out in the case 

Fig. [7J We plot the critical value m c (e) as obtained from of periodic boundary conditions, which are simpler from 

the cluster coefficient Gm (as stated before, m c is defined the mathematical point of view. 

as the value at which Gm is maximum) under adsorbing The first observation one extracts from the Monte Carlo 

boundary conditions together with the critical lines pre- simulations is the dispersion in the opinions within oth- 
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erwise well defined opinion clusters. We believe that this simulations and the master equation agree initially, but 

effect comes closer to realistic situations in which a popu- start to deviate after a time that depends on the num- 

lation splits into different groups, but the opinions of indi- ber of individuals N: the smaller the size, the earlier the 

viduals within each group are not identical to each other, deviation occurs. For example, it is possible to observe 

The detailed dynamics of those clusters depends strongly in the simulations bistability between one state with only 

on the average size of the random jumps 7. For small val- one cluster (full consensus amongst the entire population) 

ues of 7, the center of mass of each cluster performs a and another state with two clusters (lack of consensus or 

random walk through the whole opinion space with an polarization in the opinions), whereas this does not ap- 

effective diffusion coefficient that scales with 1/N. As a pear in the master equation. We attribute this difference 

consequence, opinion clusters start to collide and merge to the inherent fluctuations arising from the finite number 

in a coarsening process that leads to a single large cluster of individuals of the Monte Carlo simulations. Of course, 

at very long times. For large values of 7, the mobility of for practical applications, the number of individuals will 

clusters is reduced and several opinion clusters can appear, be always finite and hence the predictions of the Monte 

Thus, a small diffusion mobility favors opinion consensus. Carlo simulations should be more relevant than those of 

the master equation. 



An order-disorder transition to cluster formation in- 



Wc have derived a master equation for the probability 

density function P[x, t) which determines the individu- 

i 1 , .i at duced by diffusion has been characterized using the so- 

als density or distribution m the opinion space. JNumen- ■' 

, . c , r -r called cluster coefficient Gm- Large 7 or to, or small e, lead 

cal integration 01 tins equation starting from uniform mi- ' ' ' 

, . , , , r nl , . . to disordered (quasihomogeneous or non-clustered state), 

tial conditions reveals that tor small 7 only one opinion v ' 

, , , „ _ . c t r 1 The value of Gm in the disordered state strongly depends 

cluster centered around x — 0.5 is formed for aimost any o j r 

, r , , . AT r-i 1 on the parameter 7. For small 7, G a 1 reaches a maximum 

value of the parameter e, as observed m the Monte Carlo ^ ' " 

, , n 1 r u-r value close to 1. The presence of individual's opinions ad- 

snnulations. .bor large 7, a sequence 01 bifurcations can 1 ^ 

, n , . , . , , . . , sorbcd at the extremes and 1 does not disturb signifi- 

be obtained m the opinion space. 1ms is similar to what 

r , . . . , n n , , . cantly the formation of an otherwise homogeneous state, 

was found m the original Deliuant et al. model or m our J 

j ■ pr^irm m v. 1 c j +u 4 For 7 large enough, the number of opinions adsorbed by 
previous modification [13 , 14j . We have also found that ' b ° ^ J 

when diffusion is included, the asymptotic steady-state 



the extremes increases and Gm becomes more sensitive to 



, , ,. , ., ,. , , , , , , ^ , . to and saturates far away from its maximum value 1. 

probability distributions reached by Monte-Carlo simu- 

lations are not always well represented by the ones ob- We have presented a linear stability analysis that as- 

tained from the master equation dynamics starting from sumes periodic boundary conditions at the ends of the 

the same symmetric initial condition. The Monte Carlo [0, 1] interval. This analysis allows us to derive approxi- 
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